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Abstract 

We propose and study a model of scale-free growing networks that gives a degree distribution 
dominated by a power-law behavior with a model-dependent, hence tunable, exponent. The model 
represents a hybrid of the growing networks based on popularity-driven and fitness-driven prefer- 
ential attachments. As the network grows, a newly added node establishes m new links to existing 
nodes with a probability p based on popularity of the existing nodes and a probability 1—p based 
on fitness of the existing nodes. An explicit form of the degree distribution P(p, k) is derived 
within a mean field approach. For reasonably large k, P(p,k) ~ fc _7 ( p ).F(/c,p), where the function 
J- is dominated by the behavior of l/ln(/c/m) for small values of p and becomes ^-independent as 
p — ► 1, and j(p) is a model-dependent exponent. The degree distribution and the exponent 7(p) 
are found to be in good agreement with results obtained by extensive numerical simulations. 

PACS numbers: 89.75.Hc, 05.65.+b, 02.50.Cw 
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Many complex systems, including social, biological, physical, economic, and computer 
systems, can be studied by network models in which the nodes represent the constituents 
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and links or edges represent the interactions between constituents 
findings in the statistics of real-world networks reveal that classical random networks 
do not often represent the geometrical or topological structure of real-world networks 
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so-called "six degree of separation" |2Cj from one node to any arbitrary node, and the 
highly clustering feature. In particular, many networks show a power-law degree distribution 
of the form P(k) ~ /c~ 7 , with the exponent 7 taking on values between 2 to 3. This behavior 
has led to the construction of models of scale-free growing networks. Barabasi and Albert 
(BA) 21 1 proposed a model in which a new node is added in each turn and m new links are 
established with existing nodes with the probability of establishing a link being proportional 
to the number of existing links of the nodes. This preferential attachment is thus driven 
entirely by the popularity of existing nodes. Detailed numerical simulations and analytic 
analysis showed that 7 = 3 for the BA model. The idea of incorporating preferential 
attachment in a growing network has led to proposals of a considerable number of models 
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Alternativel y m odels with preferential attachment 
driven entirely by fitness have also been proposed |26l.l27|. In these models, each node carries 
a randomly assigned fitness value that gives a collective character of the node other than its 
popularity, and the probability of establishing a new link to an existing node is proportional 
to the product of fitness and the number of existing links. It was found that the degree 
distribution P(k) ~ A; -7 / ln(fc/m), with 7 = 2.255 for the fitness-driven model 26]. 

Many real-world networks such as paper citations in scientific journals, the world-wide 
web, the internet, and the collaborative networks of actors and actresses, exhibit a degree 
distribution with a network-dependent exponent that takes on values close to but below 3. It 
is therefore interesting to construct and analyze models with a tunable degree-distribution 
exponent. In the present work, we propose and study a model representing a hybrid of 
the growing network models based on popularity-driven and fitness-driven preferential at- 
tachments of new links. As the network grows, a newly added node has a probability p of 
being popularity-driven and a probability 1 — p of being fitness-driven in establishing new 
links. Thus the resulting network consists of a mixture of two types of nodes, with a frac- 
tion p establishing new links based entirely on popularity consideration. The model reflects 
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the fact that not every one, taking the nodes as agents in a population, prefers to follow 
the popular persons, but instead may prefer to establish relationships with others based 
on characters other than the popularity of the agents. Our model thus incorporates the 
inhomogeneous nature of many real-world networks in which not all the nodes are identical. 
We report results of extensive numerical simulations on the degree distribution for networks 
of size 10 7 nodes, and compare numerical results with an analytic expression derived via 
a mean field approach. For reasonably large k, the degree distribution P(p, k) follows the 
form k~ 1 ^J-'(k,p), where J-(k,p) is dominated by the behavior of l/\n(k/m) for small p 
and becomes ^-independent for p — > 1. The exponent j(p) can be extracted numerically 
and results are found to be in good agreement with that of the mean field theory. 

Our model is defined as follows. Initially a fully connected network of mo nodes is 
constructed, with mo typically of order unity. In this work, we use mo = 5. The network 
grows with one new node being added to the existing network at a time. Each newly added 
node establishes a number of m new links to existing nodes. With probability p, the new 
node establishes links by preferential attachments based on popularity of the existing nodes 
[^l, i.e., the probability that an existing node i is connected is proportional to the degree 
or the number of links that node % carries. With probability 1 — p, the new node 
establishes links by preferential attachments based on fitness of the existing nodes j^J, i.e., 
the probability that an existing node i is connected is proportional to the product ki(t)rji, 
where the fitness 77 of a node is a randomly assigned value in the interval < i] < 1 associated 
with the node when it is introduced into the network. For y =1 (p = 0), the present model 
reduces to the popularity-driven BA [2]] (fitness-driven j^y]) model. 

Detailed numerical simulations have been carried out for our model with networks of 
size up to 10 7 nodes. Each newly added node establishes m = 5 new links. Figure 1 
shows a typical degree distribution function on a log-log scale for the case of p = 0.5, 
i.e., a new node randomly chooses to follow the popularity-driven or fitness-driven rules in 
establishing new links, and for p = 0.2 (see inset). The large size of the networks used in 
this study makes the comparison with analytic results and the extraction of the exponent 
in the degree distribution easier. The data shown in Fig.l represent an average over 10 
different realizations of networks of the same size. To explore the functional form of the 
degree distribution and to extract possible exponent, we need some guidance from analytic 
treatment. 
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The model can be treated analytically by a mean field approach |1L l2t 1261] . For sufficiently 
long time, the connectivity ki(t) of the i-th node with fitness T]i evolves according to the 
following continuous time evolution equation 

dki kAt) ,„ . riikAt) . 

- =mp ^l + m(l - p) ±A±, (1) 

where the first and second terms describe the increment in ki due to popularity and fitness, 
respectively. Note that J2j kj = 2mt, with the factor of 2 coming from the undirected nature 
of the links. We assume that ki takes on the form 

/ f \ P{m,v) 

ki(t,to) = m - , (2) 



do J 

where to is the time at which the i-th node was introduced into the network. Since ki can 
at most be increased by one at each time step, it cannot grow faster than t, thus implying 
0</%,p)<l. 

The fitness rji is, in general, chosen from a distribution p(rj). The average of the sum 
J2jVjkj over P(v) can be evaluated by 

rt 



(Y.Vjkj) = J VjPiVj) J k j (t,t )dt dr] j 



For large t, since t^/t — > 0, the contribution from the term t^ r,i ' p > becomes negligible and 
we have 

(5>*i> = C(p)m*, (4) 



where 



Substituting Eqs.(2) and (4) into Eq.(l), we obtain 

^p) = l + ^Vi- (6) 

From now on, we drop the subscripts i and j for brevity. Equation (6) and Eq.(5) can be 
solved self-consistently for C(p). For p(rj) being a uniform distribution between and 1, i.e., 
p(rj) = 1 for 7] G [0, 1], C{jp) satisfies 

^ 2 ^ C(p) 



for which the integral can be performed to obtain the self-consistent equation 

1 l-p/2 A 2d-p) \ 
— — — r = 1 - e c (P ) . g) 

C(p) 1 - p V J K J 

Equation (8) can be solved numerically for C(p), for given value of p. It is found that 
C(p) G [1, 1.255] for p G [0, 1]. Note that for p = 1, ft = 1/2 (see Eq.(6)) independent of 77 
as in the BA popularity- driven mode l 12 ill . For p = 0, C(0) = 1.255 and ^(77, 0) = 77/1.255, 
as in the purely fitness-driven model 26]. 

To proceed, the cumulative probability distribution function (CDF) V v (p,ki > k) for a 
particular fitness 77 and given p can be found by noting that a degree higher than some 
value k for a node corresponds to a cutoff in time before which the node must have been 
introduced into the network, i.e., 

V v (p,ki > k) — V n {t G <t[-\ ) 

= (t) • (9) 

Note that a pre-factor of t/(rriQ + 1), which approaches unity for sufficiently long time, has 
been ignored in Eq.(9). To obtain the CDF for the whole network V(P,k), an average is 
taken over a uniform distribution of rj, thus 

V{p,k)= / - d v , (10) 



k 

with j3(r],p) given by Eq.(6). The probability distribution function (PDF) of degrees in the 
network P(p, k) is related to V(p, k) through 

Pfak) = ~P(p,k). (11) 

Using Eq.(10) and by making the substitution x = ln(/c/m)/(p/2 + (1 —p)i]/C(p)), Eq.(ll) 
gives 

^■*)=*&jC^ (12) 

where the lower limit of the integral xq = aoln(k/m) and the upper limit x\ = an \n(k/m), 
with 

a ° = p/2 + (l-p)/C(p) (13) 

and 

a x = -. (14) 
p 



We aim at getting the functional form of P(p, k) and in particular the dependence on 
k at given p. The integral in Eq.(12) can be carried out by parts, and only the "surface" 
term survives in the limit of ln(k/m) 3> 1. It follows from the fact that 1/x takes on the 
maximum value of x^ 1 for x in the interval [scoj^i], and thus the integral J^(e~ x /x 2 )dx < 
Xq 1 J^(e~ x /x)dx, with Xq 1 <C 1 for large k. Equation (12) thus gives 

m(l — p) a m{k/m) I a,\ \m J I 

Equation (15) is the main result of the mean field treatment. It gives the explicit k- 
dependence of the degree distribution function. It is worth noting that Eq.(15) gives 
the correct results in both limits of y — » and p — > 1. For purely fitness-driven model, 
P(0,k) = (k I m) -2,255 / (mln(k / m)) 26j. For purely popularity-driven model, the \n(k/m) 
term in the denominator can be shown to be cancelled by the terms in the parentheses, 
giving P(l, k) = 2m 2 k~ 3 Q- 

Results obtained from the mean field theory can be compared with numerical results. The 
solid lines in Fig.l show the degree distribution for p = 0.5 andp = 0.2 (inset) using Eq.(12). 
Excellent agreements are found between mean field and numerical results. Equation (15) 
also suggests a functional form for P(p, k). For reasonably large k, e.g., k ~ 10 2 or above, 
the second term in the parentheses is small compared to unity, especially for small values of 
p. Hence, Eq.(15) suggests that P(p,k) ~ k~ 7 ^ / \n(k / m) , where 

j(p) = 1 + a = 1 + — — — ^ , 7777 (16) 

( p /2 + (l-p)/C(p)) 

with C(p) given by Eq.(8) within mean field theory. Numerically, we fit the degree distri- 
bution function to the form fc -7 ^/ \n(k/m) and extract the exponent 7Q9) directly from 
results of numerical simulations for each value of p in steps of 0.1. However, as p — > 1, one 
should be more careful in handling Eq.(15) in numerical extraction of the exponent as the 
prefactor 1/(1 — p) diverges and the \n(k/m) term becomes unimportant due to cancellation 
effect from the terms in the parentheses. In this case, it is more convenient to work from 
the p — 1 limit and extract a functional form from Eq.(15) that is valid for p — > 1. The 
result is P(p, k) ~ (p/2 + (1 — p)/C)(k/m)~ 1 ^ /mp 2 , where 7(p) is again given by Eq.(16). 
This form is used to extract the exponent j(p) for p = 0.9 and 1.0. For p < 0.8, using either 
functional form extracts the same value of 7. Figure 2 shows the exponent ^(p) numerically 
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extracted from simulations, together with the analytic result. The two sets of results are 
found to be in good agreement. 

In summary, we proposed and studied a model in which the nodes are inhomogeneous. 
The model combines popularity-driven and fitness-driven preferential attachments in grow- 
ing networks. Extensive numerical simulations were carried out and a mean field theory was 
developed. The degree distribution function shows a predominant power-law behavior. The 
exponent takes on a model-dependent, hence tunable, value depending on the concentration 
p of nodes for which the links are established by a popularity-driven mechanism. The ex- 
ponent j(p) takes on values between 2.255 and 3, which lie within the range of values of 
7 observed in many real-world networks. Analytic expressions for the degree distribution 
function and the exponent 7(p) were derived. Results of mean field theory were found to be 
in good agreement with results obtained by numerical simulations. 
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FIGURE CAPTIONS 



Figure 1: The degree distribution function P(p, k) as a function of A; on a log-log scale for 
p = 0.5 and p = 0.2 (inset). The symbols give numerical results averaging over 10 different 
realizations of networks of size of 10 7 nodes. The lines give the analytic results within the 
mean field theory. 

Figure 2: The exponent j(p) characterizing the degree distribution function as obtained by 
fitting to numerical simulation results (symbols) for networks of size 10 7 nodes for values of 
p from p = to p = 1 in steps of 0.1, and by the mean field theory (dotted line) given by 
7 (p) = 1 + l/(p/2 + (1 -p)/C{p)) with C{p) given by Eq.(8). 
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